ARID1A regulates DNA repair through chromatin organization and its deficiency triggers DNA damage-mediated anti-tumor immune response.

AT-rich interaction domain protein 1A (ARID1A), a SWI/SNF chromatin remodeling complex subunit, is frequently mutated across various cancer entities. Loss of ARID1A leads to DNA repair defects. Here, we show that ARID1A plays epigenetic roles to promote both DNA double-strand breaks (DSBs) repair pathways, non-homologous end-joining (NHEJ) and homologous recombination (HR). ARID1A is accumulated at DSBs after DNA damage and regulates chromatin loops formation by recruiting RAD21 and CTCF to DSBs. Simultaneously, ARID1A facilitates transcription silencing at DSBs in transcriptionally active chromatin by recruiting HDAC1 and RSF1 to control the distribution of activating histone marks, chromatin accessibility, and eviction of RNAPII. ARID1A depletion resulted in enhanced accumulation of micronuclei, activation of cGAS-STING pathway, and an increased expression of immunomodulatory cytokines upon ionizing radiation. Furthermore, low ARID1A expression in cancer patients receiving radiotherapy was associated with higher infiltration of several immune cells. The high mutation rate of ARID1A in various cancer types highlights its clinical relevance as a promising biomarker that correlates with the level of immune regulatory cytokines and estimates the levels of tumor-infiltrating immune cells, which can predict the response to the combination of radio- and immunotherapy.

Dnmt3a1 regulates hippocampus-dependent memory via the downstream target Nrp1.

Epigenetic factors are well-established players in memory formation. Specifically, DNA methylation is necessary for the formation of long-term memory in multiple brain regions including the hippocampus. Despite the demonstrated role of DNA methyltransferases (Dnmts) in memory formation, it is unclear whether individual Dnmts have unique or redundant functions in long-term memory formation. Furthermore, the downstream processes controlled by Dnmts during memory consolidation have not been investigated. In this study, we demonstrated that Dnmt3a1, the predominant Dnmt in the adult brain, is required for long-term spatial object recognition and contextual fear memory. Using RNA sequencing, we identified an activity-regulated Dnmt3a1-dependent genomic program in which several genes were associated with functional and structural plasticity. Furthermore, we found that some of the identified genes are selectively dependent on Dnmt3a1, but not its isoform Dnmt3a2. Specifically, we identified Neuropilin 1 (Nrp1) as a downstream target of Dnmt3a1 and further demonstrated the involvement of Nrp1 in hippocampus-dependent memory formation. Importantly, we found that Dnmt3a1 regulates hippocampus-dependent memory via Nrp1. In contrast, Nrp1 overexpression did not rescue memory impairments triggered by reduced Dnmt3a2 levels. Taken together, our study uncovered a Dnmt3a-isoform-specific mechanism in memory formation, identified a novel regulator of memory, and further highlighted the complex and highly regulated functions of distinct epigenetic regulators in brain function.

Discovery of non-genomic drivers of YAP signaling modulating the cell plasticity in CRC tumor lines.

In normal intestines, a fetal/regenerative/revival cell state can be induced upon inflammation. This plasticity in cell fate is also one of the current topics in human colorectal cancer (CRC). To dissect the underlying mechanisms, we generated human CRC organoids with naturally selected genetic mutation profiles and exposed them to two different conditions by modulating the extracellular matrix (ECM). Among tested mutation profiles, a fetal/regenerative/revival state was induced following YAP activation via a collagen type I-enriched microenvironment. Mechanistically, YAP transcription was promoted by activating AP-1 and TEAD-dependent transcription and suppressing intestinal lineage-determining transcription via mechanotransduction. The phenotypic conversion was also involved in chemoresistance, which could be potentially resolved by targeting the underlying YAP regulatory elements, a potential target of CRC treatment.

Aberrant MNX1 expression associated with t(7;12)(q36;p13) pediatric acute myeloid leukemia induces the disease through altering histone methylation.

Certain subtypes of acute myeloid leukemia (AML) in children have inferior outcome, such as AML with translocation t(7;12)(q36;p13) leading to an MNX1::ETV6 fusion along with high expression of MNX1. We have identified the transforming event in this AML and possible ways of treatment. Retroviral expression of MNX1 was able to induce AML in mice, with similar gene expression and pathway enrichment to t(7;12) AML patient data. Importantly, this leukemia was only induced in immune incompetent mice using fetal but not adult hematopoietic stem and progenitor cells. The restriction in transforming capacity to cells from fetal liver is in alignment with t(7;12)(q36;p13) AML being mostly seen in infants. Expression of MNX1 led to increased histone 3 lysine 4 mono-, di- and trimethylation, reduction in H3K27me3, accompanied with changes in genome-wide chromatin accessibility and genome expression, likely mediated through MNX1 interaction with the methionine cycle and methyltransferases. MNX1 expression increased DNA damage, depletion of the Lin-/Sca1+/c-Kit+ population and skewing toward the myeloid lineage. These effects, together with leukemia development, were prevented by pre-treatment with the S-adenosylmethionine analog Sinefungin. In conclusion, we have shown the importance of MNX1 in development of AML with t(7;12), supporting a rationale for targeting MNX1 and downstream pathways.

Clinical impact of the genomic landscape and leukemogenic trajectories in non-intensively treated elderly acute myeloid leukemia patients.

To characterize the genomic landscape and leukemogenic pathways of older, newly diagnosed, non-intensively treated patients with AML and to study the clinical implications, comprehensive genetics analyses were performed including targeted DNA sequencing of 263 genes in 604 patients treated in a prospective Phase III clinical trial. Leukemic trajectories were delineated using oncogenetic tree modeling and hierarchical clustering, and prognostic groups were derived from multivariable Cox regression models. Clonal hematopoiesis-related genes (ASXL1, TET2, SRSF2, DNMT3A) were most frequently mutated. The oncogenetic modeling algorithm produced a tree with five branches with ASXL1, DDX41, DNMT3A, TET2, and TP53 emanating from the root suggesting leukemia-initiating events which gave rise to further subbranches with distinct subclones. Unsupervised clustering mirrored the genetic groups identified by the tree model. Multivariable analysis identified FLT3 internal tandem duplications (ITD), SRSF2, and TP53 mutations as poor prognostic factors, while DDX41 mutations exerted an exceptionally favorable effect. Subsequent backwards elimination based on the Akaike information criterion delineated three genetic risk groups: DDX41 mutations (favorable-risk), DDX41wildtype/FLT3-ITDneg/TP53wildtype (intermediate-risk), and FLT3-ITD or TP53 mutations (high-risk). Our data identified distinct trajectories of leukemia development in older AML patients and provide a basis for a clinically meaningful genetic outcome stratification for patients receiving less intensive therapies.

Heterogeneity in the gene regulatory landscape of leiomyosarcoma.

Characterizing inter-tumor heterogeneity is crucial for selecting suitable cancer therapy, as the presence of diverse molecular subgroups of patients can be associated with disease outcome or response to treatment. While cancer subtypes are often characterized by differences in gene expression, the mechanisms driving these differences are generally unknown. We set out to model the regulatory mechanisms driving sarcoma heterogeneity based on patient-specific, genome-wide gene regulatory networks. We developed a new computational framework, PORCUPINE, which combines knowledge on biological pathways with permutation-based network analysis to identify pathways that exhibit significant regulatory heterogeneity across a patient population. We applied PORCUPINE to patient-specific leiomyosarcoma networks modeled on data from The Cancer Genome Atlas and validated our results in an independent dataset from the German Cancer Research Center. PORCUPINE identified 37 heterogeneously regulated pathways, including pathways representing potential targets for treatment of subgroups of leiomyosarcoma patients, such as FGFR and CTLA4 inhibitory signaling. We validated the detected regulatory heterogeneity through analysis of networks and chromatin states in leiomyosarcoma cell lines. We showed that the heterogeneity identified with PORCUPINE is not associated with methylation profiles or clinical features, thereby suggesting an independent mechanism of patient heterogeneity driven by the complex landscape of gene regulatory interactions.

The GAPDH redox switch safeguards reductive capacity and enables survival of stressed tumour cells.

Glyceraldehyde 3-phosphate dehydrogenase (GAPDH) is known to contain an active-site cysteine residue undergoing oxidation in response to hydrogen peroxide, leading to rapid inactivation of the enzyme. Here we show that human and mouse cells expressing a GAPDH mutant lacking this redox switch retain catalytic activity but are unable to stimulate the oxidative pentose phosphate pathway and enhance their reductive capacity. Specifically, we find that anchorage-independent growth of cells and spheroids is limited by an elevation of endogenous peroxide levels and is largely dependent on a functional GAPDH redox switch. Likewise, tumour growth in vivo is limited by peroxide stress and suppressed when the GAPDH redox switch is disabled in tumour cells. The induction of additional intratumoural oxidative stress by chemo- or radiotherapy synergized with the deactivation of the GAPDH redox switch. Mice lacking the GAPDH redox switch exhibit altered fatty acid metabolism in kidney and heart, apparently in compensation for the lack of the redox switch. Together, our findings demonstrate the physiological and pathophysiological relevance of oxidative GAPDH inactivation in mammals.

TAL1 activation in T-cell acute lymphoblastic leukemia: a novel oncogenic 3' neo-enhancer.

T-cell acute lymphocytic leukemia protein 1 (TAL1) is one of the most frequently deregulated oncogenes in T-cell acute lymphoblastic leukemia (T-ALL). Its deregulation can occur through diverse cis-alterations, including SIL-TAL1 microdeletions, translocations with T-cell Receptor loci, and more recently described upstream intergenic non-coding mutations. These mutations consist of recurrent focal microinsertions that create an oncogenic neo-enhancer accompanied by activating epigenetic marks. This observation laid the groundwork for an innovative paradigm concerning the activation of proto-oncogenes via genomic alterations of non-coding intergenic regions. However, for the majority of T-ALL expressing TAL1 (TAL1+), the deregulation mechanism remains 'unresolved'. We took advantage of H3K27ac and H3K4me3 chromatin immunoprecipitation sequencing data of eight cases of T-ALL, including five TAL1+ cases. We identified a putative novel oncogenic neo-enhancer downstream of TAL1 in an unresolved monoallelic TAL1+ case. A rare but recurrent somatic heterozygous microinsertion within this region creates a de novo binding site for MYB transcription factor. Here we demonstrate that this mutation leads to increased enhancer activity, gain of active epigenetic marks, and TAL1 activation via recruitment of MYB. These results highlight the diversity of non-coding mutations that can drive oncogene activation.

Dynamic DNA methylation reveals novel cis-regulatory elements in mouse hematopoiesis.

Differentiation of hematopoietic stem and progenitor cells to terminally differentiated immune cells is accompanied by large-scale remodeling of the DNA methylation landscape. Although significant insights into the molecular mechanisms of hematopoietic tissue regeneration were derived from mouse models, profiling of DNA methylation has been hampered by high cost or low resolution using available methods. The recent development of the Infinium Mouse Methylation BeadChip (MMBC) array facilitates methylation profiling of the mouse genome at a single CpG resolution at affordable cost. We extended the RnBeads package to provide a computational framework for the analysis of MMBC data. This framework was applied to a newly generated reference map of mouse hematopoiesis encompassing nine different cell types. Analysis of dynamically regulated CpG sites showed progressive and unidirectional DNA methylation changes from hematopoietic stem and progenitor cells to differentiated hematopoietic cells and allowed the identification of lineage- and cell type-specific DNA methylation programs. Comparison with previously published catalogs of cis-regulatory elements (CREs) revealed 12,856 novel putative CREs that were dynamically regulated by DNA methylation (mdCREs). These mdCREs were predominantly associated with patterns of cell type-specific DNA hypomethylation and could be identified as epigenetic control regions regulating the expression of key hematopoietic genes during differentiation. In summary, we established an analysis pipeline for MMBC data sets and provide a DNA methylation atlas of mouse hematopoiesis. This resource allowed us to identify novel putative CREs involved in hematopoiesis and will serve as a platform to study epigenetic regulation of normal and malignant hematopoiesis.

Global DNA Methylation Analysis of Cancer-Associated Fibroblasts Reveals Extensive Epigenetic Rewiring Linked with RUNX1 Upregulation in Breast Cancer Stroma.

Cancer cells recruit and rewire normal fibroblasts in their microenvironment to become protumorigenic cancer-associated fibroblasts (CAF). These CAFs are genomically stable, yet their transcriptional programs are distinct from those of their normal counterparts. Transcriptional regulation plays a major role in this reprogramming, but the extent to which epigenetic modifications of DNA also contribute to the rewiring of CAF transcription is not clear. Here we address this question by dissecting the epigenetic landscape of breast CAFs. Applying tagmentation-based whole-genome bisulfite sequencing in a mouse model of breast cancer, we found that fibroblasts undergo massive DNA methylation changes as they transition into CAFs. Transcriptional and epigenetic analyses revealed RUNX1 as a potential mediator of this process and identified a RUNX1-dependent stromal gene signature. Coculture and mouse models showed that both RUNX1 and its stromal signature are induced as normal fibroblasts transition into CAFs. In breast cancer patients, RUNX1 was upregulated in CAFs, and expression of the RUNX1 signature was associated with poor disease outcome, highlighting the relevance of these findings to human disease. This work presents a comprehensive genome-wide map of DNA methylation in CAFs and reveals a previously unknown facet of the dynamic plasticity of the stroma. SIGNIFICANCE: The first genome-wide map of DNA methylation in breast cancer-associated fibroblasts unravels a previously unknown facet of the dynamic plasticity of the stroma, with far-reaching therapeutic implications.

Genetic and Methylation Analysis of CTNNB1 in Benign and Malignant Melanocytic Lesions.

Melanocytic neoplasms have been genetically characterized in detail during the last decade. Recurrent CTNNB1 exon 3 mutations have been recognized in the distinct group of melanocytic tumors showing deep penetrating nevus-like morphology. In addition, they have been identified in 1-2% of advanced melanoma. Performing a detailed genetic analysis of difficult-to-classify nevi and melanomas with CTNNB1 mutations, we found that benign tumors (nevi) show characteristic morphological, genetic and epigenetic traits, which distinguish them from other nevi and melanoma. Malignant CTNNB1-mutant tumors (melanomas) demonstrated a different genetic profile, instead grouping clearly with other non-CTNNB1 melanomas in methylation assays. To further evaluate the role of CTNNB1 mutations in melanoma, we assessed a large cohort of clinically sequenced melanomas, identifying 38 tumors with CTNNB1 exon 3 mutations, including recurrent S45 (n = 13, 34%), G34 (n = 5, 13%), and S27 (n = 5, 13%) mutations. Locations and histological subtype of CTNNB1-mutated melanoma varied; none were reported as showing deep penetrating nevus-like morphology. The most frequent concurrent activating mutations were BRAF V600 (n = 21, 55%) and NRAS Q61 (n = 13, 34%). In our cohort, four of seven (58%) and one of nine (11%) patients treated with targeted therapy (BRAF and MEK Inhibitors) or immune-checkpoint therapy, respectively, showed disease control (partial response or stable disease). In summary, CTNNB1 mutations are associated with a unique melanocytic tumor type in benign tumors (nevi), which can be applied in a diagnostic setting. In advanced disease, no clear characteristics distinguishing CTNNB1-mutant from other melanomas were observed; however, studies of larger, optimally prospective, cohorts are warranted.

Club cells employ regeneration mechanisms during lung tumorigenesis.

The high plasticity of lung epithelial cells, has for many years, confounded the correct identification of the cell-of-origin of lung adenocarcinoma (LUAD), one of the deadliest malignancies worldwide. Here, we employ lineage-tracing mouse models to investigate the cell of origin of Eml4-Alk LUAD, and show that Club and Alveolar type 2 (AT2) cells give rise to tumours. We focus on Club cell originated tumours and find that Club cells experience an epigenetic switch by which they lose their lineage fidelity and gain an AT2-like phenotype after oncogenic transformation. Single-cell transcriptomic analyses identified two trajectories of Club cell evolution which are similar to the ones used during lung regeneration, suggesting that lung epithelial cells leverage on their plasticity and intrinsic regeneration mechanisms to give rise to a tumour. Together, this study highlights the role of Club cells in LUAD initiation, identifies the mechanism of Club cell lineage infidelity, confirms the presence of these features in human tumours, and unveils key mechanisms conferring LUAD heterogeneity.

Systematic evaluation of cell-type deconvolution pipelines for sequencing-based bulk DNA methylomes.

DNA methylation analysis by sequencing is becoming increasingly popular, yielding methylomes at single-base pair and single-molecule resolution. It has tremendous potential for cell-type heterogeneity analysis using intrinsic read-level information. Although diverse deconvolution methods were developed to infer cell-type composition based on bulk sequencing-based methylomes, systematic evaluation has not been performed yet. Here, we thoroughly benchmark six previously published methods: Bayesian epiallele detection, DXM, PRISM, csmFinder+coMethy, ClubCpG and MethylPurify, together with two array-based methods, MeDeCom and Houseman, as a comparison group. Sequencing-based deconvolution methods consist of two main steps, informative region selection and cell-type composition estimation, thus each was individually assessed. With this elaborate evaluation, we aimed to establish which method achieves the highest performance in different scenarios of synthetic bulk samples. We found that cell-type deconvolution performance is influenced by different factors depending on the number of cell types within the mixture. Finally, we propose a best-practice deconvolution strategy for sequencing data and point out limitations that need to be handled. Array-based methods-both reference-based and reference-free-generally outperformed sequencing-based methods, despite the absence of read-level information. This implies that the current sequencing-based methods still struggle with correctly identifying cell-type-specific signals and eliminating confounding methylation patterns, which needs to be handled in future studies.

Genetic and methylation profiles distinguish benign, malignant and spitzoid melanocytic tumors.

Accurate classification of melanocytic tumors is important for prognostic evaluation, treatment and follow-up protocols of patients. The majority of melanocytic proliferations can be classified solely based on clinical and pathological criteria, however in select cases a definitive diagnostic assessment remains challenging and additional diagnostic biomarkers would be advantageous. We analyzed melanomas, nevi, Spitz nevi and atypical spitzoid tumors using parallel sequencing (exons of 611 genes and 507 gene translocation analysis) and methylation arrays (850k Illumina EPIC). By combining detailed genetic and epigenetic analysis with reference-based and reference-free DNA methylome deconvolution we compared Spitz nevi to nevi and melanoma and assessed the potential for these methods in classifying challenging spitzoid tumors. Results were correlated with clinical and histologic features. Spitz nevi were found to cluster independently of nevi and melanoma and demonstrated a different mutation profile. Multiple copy number alterations and TERT promoter mutations were identified only in melanomas. Genome-wide methylation in Spitz nevi was comparable to benign nevi while the Leukocytes UnMethylation for Purity (LUMP) algorithm in Spitz nevi was comparable to melanoma. Histologically difficult to classify Spitz tumor cases were assessed which, based on methylation arrays, clustered between Spitz nevi and melanoma and in terms of genetic profile or copy number variations demonstrated worrisome features suggesting a malignant neoplasm. Comprehensive sequencing and methylation analysis verify Spitz nevi as an independent melanocytic entity distinct from both nevi and melanoma. Combined genetic and methylation assays can offer additional insights in diagnosing difficult to classify Spitzoid tumors.

A leukemia-protective germline variant mediates chromatin module formation via transcription factor nucleation.

Non-coding variants coordinate transcription factor (TF) binding and chromatin mark enrichment changes over regions spanning >100 kb. These molecularly coordinated regions are named "variable chromatin modules" (VCMs), providing a conceptual framework of how regulatory variation might shape complex traits. To better understand the molecular mechanisms underlying VCM formation, here, we mechanistically dissect a VCM-modulating noncoding variant that is associated with reduced chronic lymphocytic leukemia (CLL) predisposition and disease progression. This common, germline variant constitutes a 5-bp indel that controls the activity of an AXIN2 gene-linked VCM by creating a MEF2 binding site, which, upon binding, activates a super-enhancer-like regulatory element. This triggers a large change in TF binding activity and chromatin state at an enhancer cluster spanning >150 kb, coinciding with subtle, long-range chromatin compaction and robust AXIN2 up-regulation. Our results support a model in which the indel acts as an AXIN2 VCM-activating TF nucleation event, which modulates CLL pathology.

Epigenomic technologies for precision oncology.

Epigenetic patterns in a cell control the expression of genes and consequently determine the phenotype of a cell. Cancer cells possess altered epigenomes which include aberrant patterns of DNA methylation, histone tail modifications, nucleosome positioning and of the three-dimensional chromatin organization within a nucleus. These altered epigenetic patterns are potential useful biomarkers to detect cancer cells and to classify tumor types. In addition, the cancer epigenome dictates the response of a cancer cell to therapeutic intervention and, therefore its knowledge, will allow to predict response to different therapeutic approaches. Here we review the current state-of-the-art technologies that have been developed to decipher epigenetic patterns on the genomic level and discuss how these methods are potentially useful for precision oncology.

Histone H3K27 demethylase KDM6A is an epigenetic gatekeeper of mTORC1 signalling in cancer.

OBJECTIVE: Large-scale genome sequencing efforts of human tumours identified epigenetic modifiers as one of the most frequently mutated gene class in human cancer. However, how these mutations drive tumour development and tumour progression are largely unknown. Here, we investigated the function of the histone demethylase KDM6A in gastrointestinal cancers, such as liver cancer and pancreatic cancer. DESIGN: Genetic alterations as well as expression analyses of KDM6A were performed in patients with liver cancer. Genetic mouse models of liver and pancreatic cancer coupled with Kdm6a-deficiency were investigated, transcriptomic and epigenetic profiling was performed, and in vivo and in vitro drug treatments were conducted. RESULTS: KDM6A expression was lost in 30% of patients with liver cancer. Kdm6a deletion significantly accelerated tumour development in murine liver and pancreatic cancer models. Kdm6a-deficient tumours showed hyperactivation of mTORC1 signalling, whereas endogenous Kdm6a re-expression by inducible RNA-interference in established Kdm6a-deficient tumours diminished mTORC1 activity resulting in attenuated tumour progression. Genome-wide transcriptional and epigenetic profiling revealed direct binding of Kdm6a to crucial negative regulators of mTORC1, such as Deptor, and subsequent transcriptional activation by epigenetic remodelling. Moreover, in vitro and in vivo genetic epistasis experiments illustrated a crucial function of Deptor and mTORC1 in Kdm6a-dependent tumour suppression. Importantly, KDM6A expression in human tumours correlates with mTORC1 activity and KDM6A-deficient tumours exhibit increased sensitivity to mTORC1 inhibition. CONCLUSION: KDM6A is an important tumour suppressor in gastrointestinal cancers and acts as an epigenetic toggle for mTORC1 signalling. Patients with KDM6A-deficient tumours could benefit of targeted therapy focusing on mTORC1 inhibition.

Epigenetic reprogramming of airway macrophages promotes polarization and inflammation in muco-obstructive lung disease.

Lung diseases, such as cystic fibrosis and COPD, are characterized by mucus obstruction and chronic airway inflammation, but their mechanistic link remains poorly understood. Here, we focus on the function of the mucostatic airway microenvironment on epigenetic reprogramming of airway macrophages (AM) and resulting transcriptomic and phenotypical changes. Using a mouse model of muco-obstructive lung disease (Scnn1b-transgenic), we identify epigenetically controlled, differentially regulated pathways and transcription factors involved in inflammatory responses and macrophage polarization. Functionally, AMs from Scnn1b-transgenic mice have reduced efferocytosis and phagocytosis, and excessive inflammatory responses upon lipopolysaccharide challenge, mediated through enhanced Irf1 function and expression. Ex vivo stimulation of wild-type AMs with native mucus impairs efferocytosis and phagocytosis capacities. In addition, mucus induces gene expression changes, comparable with those observed in AMs from Scnn1b-transgenic mice. Our data show that mucostasis induces epigenetic reprogramming of AMs, leading to changes favoring tissue damage and disease progression. Targeting these altered AMs may support therapeutic approaches in patients with muco-obstructive lung diseases.